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H 1 
MTHSDTAN : Double Precision Floating Point Tangen 16-SEP-1984 01: 4 \VAX/VMS Macro Vv04-00 Pa 
MOTs . ven Sepa obs $1:55.8¢ EATANTT eaeSeruGtancmar:1 P29 
8° 1 «TITLE MTHSDTAN ; Double Precision Floating Point Tangent routine 
0 ¢ (DTAN, DTAN 
$3 € IDENT /1-013/ File: MTHDTAN.MAR €EDIT:SBL1013 
44 5 FL ASAE SRR eo 2 EOL TLE EIEN | Nk 
7* * 
$009 7 3* COPYRIGHT (c) 1978, 1980, 1982, 1984 BY * 
44 8 ;* DIGITAL EQUIPMENT CORPORATION, MAYNARD, MASSACHUSETTS. * 
44 1 :* ALL RIGHTS RESERVED.  @ 
° ® 
0000 11 ;* THIS SOFTWARE IS FURMTSHED UNDER A LICENSE AND _ BE USED AND copreD * 
464 \¢ 3* ONLY IN eee A THE TERMS OF ate ICENSE AND WITH THE * 
0 15 ;* INCLUSION OF THE ABOV BaP resettes NOTICE. THIS FTWARE OR ANY OTHER * 
0000 14 ;* COPIES THEREOF _ hee PROVIDED OR OTHERWIS OMADE AVAILABLE TO ANY * 
0000 15 ;* OTHER PERSON. ITLE Fo AND OWNERSHIP OF THE SOFTWARE IS WHEREBY * 
4 1g :* TRANSFERRED. * 
° ® 
4 18 ;* THE INFORMATION IN THIS SOFTWARE IS SUBJECT TO CHANGE WITHOUT NOTICE * 
000 19 ;* AND SHOULD NOT BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT * 
464 :* CORPORATION. * 
3 ® 
0000 ¢ 3* DIGITAL ASSUMES NO RESPONSIBILITY FOR THE USE OR RELIABILITY OF ITS * 
8000 ? :* SOFTWARE ON EQUIPMENT WHICH IS NOT SUPPLIED BY DIGITAL. ® 
3 ® 
0000 $3 3* * 
4464 37 i can sein nena nae ine aciealaanninteiaiaaiee anette 
0000 B: 
0000 $3 5 FACILITY: MATH LIBRARY 
0000 0; 
4464 31; ABSTRACT: 
0000 ‘ : MTHSDTAN is a function which returns the double precision floating point 
0000 4 ; tangent of its double precision dag | point region argument. The call is 
44 2 3 standard call-by-reference. It JSB to ATHSDTAN_R 
0000 37 : MTHSDTAND is a function which returns the double precision floating point 
0000 3 ; tangent of its double precision pleet tag poles — argument. The call is 
0000 9 ; standard call-by-reference. It JSB to MTHSDTAND_R 
0000 40 ; 
0000 41 ;-- 
0000 4g ; 
rs ; VERSION: 01 
45 : HISTORY: 
46 3 AUTHOR: 
47 ; 
48 ; 
49 ; 
50 ; 
if 
a3 
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Be 91:99:84 
Detailed Current Edit History ~SEP-1984 11:22:54 CMTHRTL. SRC JMTHDTAN.MAR; 1 


-SBTTL HISTORY ; Detailed Current Edit History 


ALGORITHMIC DIFFERENCES FROM FP-11/C ROUTINE: none 


The result is reserved operand when DCOS(X) = 0, instead 
of largest or smallest representable floating number. 


Edit History for Version 01 of MTHSDTAN 


SOoOQoooooococooooooooooooooooQooooooooooo 
WON AOULS WN | 9 OD NAU EWN OOO NAME WN ( OOONOUE 
Tr Tee TOTO rererererererererererererererererererererererere rere 


090909090909 09 0909 09 SIN NINN NINO PAA AAOAAOO 


Error pend tog mechanism changed. Instead of having 
MTHSFLG_JACKET at the entrance, MTHS$SERROR_CONT is stored o 

+44 of the stack frame so when error happened in MTHSDSIN or ATNSDCOS 
© pessage will be hided ever,» but will get signalled in MTHSDTAN. 

MIHSSERRO changed to MTHS$$S 

MTHSSERROR_CONT changed to MTHSSSIGNAL_ CON 

MTH$_... changed to ATH 

Changed arrer handling Mechanism. Put error result in RO:R1 before 

callin ng MTHS me ber in order to oes user modify error result. 

Update version number and seoxe’ ht bot ice. JBS 16-NOV-78 

anged MTH__FLOOVEMAT to MTHSK_FLOOVEMAT. JBS 06-DEC-78 

Removed S$SRADEF macro - not needed. JBS 16-DEC-78 

Add to the PSECT Nk fe S 22-DEC-7 

Declare externals. SBL 1 

Rearrange code to handle excegt long correctly. SBL 17-May-1979 

Add JSB — point. JBS 16-AU 7 

Signal excep toes correctly. UBS 16-AUG-1979 

Correct a typo in edit JBS 17-AUG-1979 

Make CALL entry use JSB routine. Correct error handling. 

Do COS before SIN. SBL 31-O0ct-197 

Added ee ree RnR? points. RNH ge? 

Modif sles ge R7 and MTHSDTAND_R7 to JSB to MTHSDSINCOS_R7 and 

MTHSDSINCOSD_R7 instead of MTHSDSIN_R7, MTHSDCOS_R7, and 

MTHSDSIND_R7, MTHSDCOSD_R7. RNH eh AUG-81. 

- Use general mode addressing. SBL 30-Nov-1981 


a A se ‘eyo Floating Point” Tanger 16-SEP =1384 91: $s: 3 ies Macro V04-00 Page (3) 


Declarative Part of Mo -SEP MTHRTL.SRCJMTHDTAN.MAR; 1 
-SBTTL DECLARATIONS ; Declarative Part of Module 


$v 
Ge 


ION 


; INCLUDE FILES: 


EXTERNAL SYMBOLS: 


-DSABL GBL 
-EXTRN MTMSDSINCOS_R7 
sEXTRN MTHSDSINCOSD_R7 
.EXTRN MTHSK_FLOOVEMAT 
-EXTRN MTHSSSIGNAL 
.EXTRN MTHSK_FLOUNDMAT 
.EXTRN MTHSSJACKET_TST 
: EQUATED SYMBOLS: none 
: MACROS: 
SSFDEF 3; Define SF (Stack Frame) symbols 


PSECT DECLARATIONS: 


-PSECT _MTHSCODE PIC,SHR, LONG, EXE ,NOWRT 
3 program section for math routines 


0000 


OWN STORAGE: none 


CONSTANTS: 
00000004 X=4 :Position of argument from AP. 
D_SMALLEST_DEG: 
SEEDS 36S -LONG  “X2EE04365 
FBEDSIE [LONG  “*XOFBED31E 


SOOCSCOOOCSOOCSOSOOSCOOSOSOOOSOOOOOOOSOOOOOOOSOSOOOSOSSSoO -8 


SOOCSCOOCSOSOSOSOOCOOSOOOSOSOSOSOSOSOSSOOSOOOSOOSOSOSOOSSSSO $Y 
DWLOOCOCCOCOCOOSCOCOCOCOSOOOCCOO COO OOOOOOOOOOOOOoOoOO 


SOONAUES WN 9S OONA UE WIN O ODNA UNE WW S$ O OONOU Ewin 


WiRoROnoNonononononons 2 3 2 3 3 2 = 3 | BODO OOOO O0O0O0000 


a a ts a tt 


SOoOCCSOOOOOOOSOOOSOSOSOOOOOOOooOOD 


6D O00000000'GF 9E MOVAB G*MTHSSJACKET_HND, (FP) 
3; set handler address to jacket 


3; handler 
171 ax(AP), RO ; RO/R1 = argument 
Vie BSBB MTHSDTAN_R7 ; Get the tangent 
17 RET ; With result in RO-R1 
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MTHSDTAN ; Double Precision Floating Point Tangen 16-SEP-1984 01:22:34 VAX/VMS Ma v04-00 P 4 
vast MTHSDTAN = Standard Double Precision ~$6P 138% 9:88:82 EMTHRTL. SREINTHD TAN MARS 1 _— (4) 
} ¢ -SBTTL MTHSDTAN = Standard Double Precision Floating DTAN 
136 4s 
! § ; FUNCTIONAL DESCRIPTION: 
! 8 ; DTAN = double precision floating point function 
0 140 : For algorithm, see MTHSDTAN_R7. 
146 ; CALLING SEQUENCE: 
144 ; DTAN.wd.v = MTHSDTAN(x.rd.r) 
090 146 > INPUT PARAMETERS: 
44 138 3 X.ed.r Address of value of angle in radians. 
2098 130 : IMPLICIT INPUTS: none 
0008 13¢ | OUTPUT PARAMETERS: 
44 136 ; VALUE: double precision floating tangent of the argument 
0008 136 ; IMPLICIT OUTPUTS: none 
9008 138 : COMPLETION CODES: none 
9008 160 3 SIDE EFFECTS: 
$008 16¢ ; See MTHSDTAN_R7 
Bp08 163 ; 
008 164 ;--- 
She 
40FC 008 167 «ENTRY MTHSDTAN, “M<IV, R2, R3, R4, RS, Rb, R7> 
A 168 : standard call-by-reference entry 
A 169 ; disable DV (and FU), enable IV 
0 A 170 MTHSFLAG_JACKET 
QO00A 
8 1 
1 
1 
0 
5 
0017 
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~SBTTL MTHSDTAN_R7 = JSB entry point 


oo 


cro V04-00 Page ; 
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ac 
ee 
a] 
NO 


=z: 
~ 
xo 
ero 
oc 


"FUNCTIONAL DESCRIPTION: 
DTAN = JSB entry point 
Te Compute ste ope: 
° ute DSIN and DCOS. 
° £F OSS Tt oer es ve have an error. 
CALLING SEQUENCE: 
H83” REARDAN a 
INPUT PARAMETERS: 
RO / R1 contains X 
IMPLICIT INPUTS: 
NONE 
OUTPUT PARAMETERS: 
The result is the tangent of X, in double precision. 
IMPLICIT OUTPUTS: 
NONE 


SIDE EFFECTS: 


Se Ge Ge Ge Se Se Ge Ge Ge Ge Ge Se Ge Ge Ge Ge Ge Ge Ge Se Ge Ge Ge Ge Se Sse Ge Se ee ae 


MTHSDTAN_R7:: 
00000000'°GF 1 JSB G*MTHSDSINCOS_R7 : Compute DSIN(X) and DCOS(X) 
5 7 TSTD : Is DCOS(X) zero ? 
6 1 BEQL COSZER : if yes, then 7 to common error path 
50 52 66 DIvD2 R2, RO 3; Compute DSIN cos 
05 RSB 3; Return. 


$i 
5 
5 
5 
8 
b: 
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ROPOPONON) 2 2 2 9 


QAO&WYOM 0900 09 09 C9 Cd Cd Cd Cd OD Cd Cd C9 CD Cd Cd Cd OD Cd Od Cd Cd Cd Cd CD CD. 0D. C9 C9 CD 09 0D CD CD CO CDOD 0000 
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~SBTTL MTHSDTAND = Standard Double Precision Floating DTAND 


{FUNCTIONAL DESCRIPTION: 
: DTAND = double precision floating point function 
For algorithm, see MTHSDTAND_R7. 
CALLING SEQUENCE: 
DTAND.wd.v = MTHSDTAND(X.rd.r) 
INPUT PARAMETERS: 
X.rd.r Address of value of angle in degrees. 
IMPLICIT INPUTS: none 
OUTPUT PARAMETERS: 
VALUE: double precision floating tangent of the argument 
IMPLICIT OUTPUTS: none 
COMPLETION CODES: none 
SIDE EFFECTS: 
See MTHSDTAND_R7 


40FC ENTRY MTHSDTAND, “M<IV, R2, R3, R4, RS, R6, R7> 
3 standard call-by-reference entry 


3; disable DV (and FU), enable IV 


UA ATU ES BS BS SS BB EB BWW WII NII OOPONONONINNND 
Wr 


DWONAUES WN O OONAUE WN 9 OD NOAUE WN O OONOUS 


rorore 


MTHSFLAG_JACKET 


6D Q0000000'GF 9 MOVAB G*MTHSSJACKET_HND, (FP) 


set handler address to jacket 


SOSSSCSOSOOSOSOOOOSOOOSOSOSOOSOSOOOSOSOOOSOOSOOOSOSOOOOOOOOOOOOOO 


SOSOSSCOCOCOOCOOOC OOOO SOSOSOSOSOOOSOOSOO OOOO SOOoOOSOOS: 


IEP IPINIPOPINININYININININININPININININIPININIPINIPUDININIPIPIPIPIPOPIPUPIPUPUPIPUPIPUNPININNY 


WENT “FIFI “79 CDOD GOOD OD O* OS OS OS OS OS ODP APPA AAA AA AA AAA AAA A AAAS 


3; handler 
50 04 ac 70 MOVD aXx(AP), RO ; RO/R1 = argument 
10 61 BSBB  MTHSDTAND_R7 : Get the tangent 
04 RET 3 With result in RO-R1 


6 
(6) 
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MTHYDTAN ; Double Precision Floating Point Ta 16-SEP-1984 01:22:34 VAX/VMS Macro V04-00 Page 7 
1-018 ATHSDTAND_R? - JSB entry paint sees o7$EE= 1 3be 97:98:82 MTHRTL.SRCIMTHDTAN.MAR; 1 we 
be : -SBTTL MTHSDTAND_R7 = JSB entry point 
66 +4 
? 3 FUNCTIONAL DESCRIPTION: 
8 ; DTAND = JSB entry point 
6 6 1: Algorithmic steps: 
8 6 1 3; 1. Make sure t st she absolute value of the argument is greater 
0036 73; than 180/pi*2*128 
036 74; ¢: compue DSIND and DCOSD. ' 
$ ° 75; 3. If DCOSD is zero, we have an error. 
036 i: ; 4. Return DSIND / ptosp. 
0 6 28 3 CALLING SEQUENCE: 
0036 80; MOVD argument, RO 
0036 81 : JSB MTASDTAND_R7 
0056 es : INPUT PARAMETERS: 
bose A : RO / R1 contains x 
0036 Ha : IMPLICIT INPUTS: 
0036 88 ; 
oe ym 
0056 4 : OUTPUT PARAMETERS: 
base $8 ; The result is the tangent of x, in double precision. 
0036 295 : IMPLICIT OUTPUTS: 
4 6 96 ; 
0036 | 598 : sete 
S036 99 : SIDE EFFECTS: 
0036 00 ; 
0036 3} 3 NONE 
036 8 MTHSDTAND_R7:: 
52 50 76 0036 04 mMOvD RO, Re 3 Save_argument 
38 8008 8F AA 0039 05 BICW #*x8000, R ; R2/R3 = ‘argument: 
5 0300 8F BI o6 CMP #*x0300, R ; Compare iarg: with 2**-121 
Od 15 ‘3 8 BLEQ $ 3; No possible underflow, compute DCOSD 
52 BBAF 71 4 8 CMPD D_SMALLEST_DEG, R2 :; Possible underflow, use better check 
97 v 49 09 BLSS 20$ : No possible underflow, compute DCOSD 
¢ 7 48 \9 TSTD R2 ; Check for argument = 
1 12 4D 1 BNEQ UNFL 3 3f arg not 0 go to underflow logic 
50 43 4F 19 CLRD = RO. : Load RO/R1 with zero 
0 2} 1? 208 RSB 3; Return with value = 0 
"GF 1 15 : JSB G*MTHSDSINCOSD_R7 s ¢ te DSIND(X), and DCOSD(X) 
nesneen 2 8 33 16 TSTD e 3 Is bCOSD(x) zero ? 
2 1 5A 1 BEQL COSZER ; If yes, then go to common error path 
50 66 ¢ H DivD2 R2, RO 3; Compute DSIND / DCOSD 
2h 4 RSB 3; Return. 
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00000000 ' GF 
52 


OD 52 
7E 


00000000 ' GF 


7E 
00000000 ' GF 


04 
04 


So 
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02 


OF 
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Foo ISB entry point Seer obe 125584 EATANTS SeeserMotageman:1 29° (By 
33 
¢ : COMMON ERROR CODE 
8 3; Underflow; if user has FU set, signal error. Always return 0.0 
0 UNFL: 
1 MOVPSL 3 ; R2 = user's or jacket routine's PSL 
¢ CALLS #0, G*MTHSSJACKET_TST ; RO = TRUE if JSB from jacket routine 
BLBC ; branch if user d iy 
4 MOVZWL SFS$W_SAVE_PSW(FP), R2  ; get user PSL saved by CA LL 
5 10$: CLRL RO ; RO = toy LIBSSIGNAL will save in 
§ ; CHFSL_M RO/R1 so any handler 9 fixup 
BBC #6, R2, 20$ 3 has ve gts Wleating under f Low? 
38 PUSHL (SP) 3; yes, return PC from special routine 
39 MOVZBL #MTHSK_FLOUNDMAT, -(SP) ; trap code for hardware floaging under f Low 
40 3; convert to MTHS_FLOUNDMAT (32-bit VAX-11 
41 3 exception code)” 
348 CALLS #2, G*MTHSSSIGNAL 3 signal (condition, PC) 
43 20$: RSB 3 return 
$e 
45 i+ 
o9 3. ; Come here if COS value is zero. 
48 coszeR: 
49 PUSHL (SP ; Push ‘“‘caller’’ PC 
50 ret paTusK FLOOVEMAT, -(SP) ; Condition value 
51 #15, ; RO/R1 = reserved operand 
26 CALLS #2, G°MTHS$SIGNAL ; Signal an error 
27 RSB 3; Return to caller 
355 - END 


MTHSOTA 
Symbol table 


COSZER 4 4 0084 R 

D_SMALLEST a5 0000000 R 
MTHSSJACKET_H eeeeeeee XX 
MINSSOACKET= aS eeeeeere =X 
MTHSSSIGNAL~ eeereere x 
MTHSDSINCOSD_R7 ‘eeeeeeee xX 
eereeeee xX 


sa 4} R7 
MTHSDTA 00000008 RG 
000000 
000000 


THSDTAN_R7 
MTHSK_FLOOVEMAT teeeeeee 
MTHSK_FLOUNDMAT tteeeree : 
SFSW_SAVE_PSW = 00000004 


Ibs working is cae was 900 page 
1 byies 

There were 
415 source Snee were read in Pass 1, 


Macro Library name 


_$255$DU >8:CSYSLIBISTARLET.MLB;2 


producin 
9 pages of virtual memory were used to define 


Coes meron ana sao ea acon + 


00000060 R 02 
¥ = 00000004 
PSECT name Allocation 
Ss 00000000 0.) 
$ 00000000 0.) 
MTHSCODE 00000096 <¢ 150.) 
Phase Page faults CPU Time 
Initialization 33 00:00:00.09 
pennene processing 112 00:00:00.70 
s 1 121 00:00:01.46 
ans table sort 0 BF 288: 80-82 
Pass 2 74 00:00:00.90 
Symbol table output 3 B98 08 Be 
Psect synopsis output 2 00:00:00.0 
Cross-reference output ! B82 09 209-99 
Assembler run totals 34 00:00:03.2 


eam mn we oe enema eee ewe t 


dD 2 
3; Double Precision Floating Point Tangen 16-SEP-1984 01:22:34 
: ° we 6- mit 94:38:34 


Seer ere esse w + 


00:00:00.84 
BO: 00: 28- 99 


ooooooo 
ooocececccs 
r=t~t=t=t=t—t-t— 
ooosceccs 


pages) of virtual enett were used to buffer the intermediate c 

0 pages of symbol table space allocated to hold 44 non-local and 5° a symbols. 
16 object records in Pass 2 

macros. 
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MTHRTL. SR RCIM MTHDTAN.MAR; 1 (8 


| 


PSECT No. Attributes 
00 ¢« O. ; NOPIC USR CON ABS LCL 

( 1.) NOPIC USR CON ABS LCL NOSHR EXE RD WRT NOVEC BYTE 
02 ¢ 2. } PIC USR CON REL LCL SHR EXE RD NOWRT NOVEC LONG 


| 
| 
| 
| 
| 
NOSHR NOEXE NORD NOWRT NOVEC BYTE 
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»Aibeies ed ae ae hs 91:88:84 MTHRTL. SRCIMTHDTAN.MAR: 1 (8) 


th fie 
VAX-11 Macro Run Statistics 
88 GETS were required to define 4 macros. 


There were no errors, warnings or information messages. 
MACRO/ENABLE=SUPPRESSION/D1 SABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:MTHDTAN/OBJ=OBJ$:MTHDTAN MSRC$:MTHJACKET/UPDATE=(ENH$:MTHJACKET) +MSRCS: 


NT CORPORATION 
D PROPRIETARY 


AH-BT13A-SE 
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